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Abstract 

Excitation spectra and neutron single particle configurations of '^'^Mg and ^^Ne are investigated 
by using antisymmetrized molecular dynamics combined with generator coordinate method. It is 
shown that both nuclei have strongly deformed 3/2^ ground state with a 3p2/i configuration. The 
excitation spectra are qualitatively understood in terms of the Nilsson model and the calculation 
has shown the coexistence of different intruder configurations within small excitation energy. The 
calculated one neutron separation energy of ^^Ne is rather small {Sn = 250 keV) and implies a 
p-wave one neutron halo with a strongly deformed core. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



The breaking of neutron magic number N — 20 that was firstly pointed out from the 
observation of the anomalous binding energy and spin of ^^Na [1, 2] has been one of the 
great interests in nuclear physics. Based on the theoretical calculations, neighboring nuclei 
around ^^Na was named "island of inversion" [3, 4] because they are dominated by the 
intruder configuration in which neutrons are promoted into p/-shell across = 20 shell gap 
due to strong deformation [5, 6]. Since then, numerous experimental and theoretical studies 
have been devoted to the island. One of the recent experimental finding is that there are 
coexistence of "spherical and deformed shape" or "normal and intruder configurations" in 
many cases [7-10] . It is expected that the exploration of the coexistence will bring us further 
understanding of the shell structure in neutron rich N 20 region. 

The spectroscopy of odd- neutron nuclei is of importance and interest, because their spin 
and parity are related to the last neutron particle's or hole's orbit. For example, we have 
predicted the coexistence of many-particle and many-hole configurations at small excitation 
energy in ^^Mg {N — 19) [11], and some of these excited states are identified by the the 
proton knock-out reaction from ^^Al [12]. The spin-parity of the excited states are associated 
with the neutron orbits [N,n„l„j,] = [2,0,2,3/2], [3,3,0,1/2], [2,0,0,1/2] and [3,2,1,3/2] in 
terms of the Nilsson model. 

In this paper, we extend our study to A'" = 21 system (^^Mg and ^^Ne). Both of them 
do not have the definite spin-parity assignment for the ground states and information of 
their excitation spectra is deficient. In the case of ^^Mg, there are two major spin-parity 
assignments of the ground state, 3/2+ and 3/2 ~. Prom the observation of the /3 decays of 
^^Na [13] and ^^Mg [14], 3/2+ assignment is suggested. On the other hand, the measurement 
of the magnetic moment [15] and one neutron knock-out reaction from '^"^Mg [16] suggest 
3/2~. ^^Ne is more interesting. Its large reaction cross section [17] and one neutron Coulomb 
breakup cross section [18] are reported. From the analysis of the one neutron breakup, the 
assignments of 1/2+ or more likely 3/2~ are suggested. In both assignments, observed large 
Coulomb breakup cross section strongly suggests one neutron-halo structure of ^^Ne that is 
firstly observed in the island. 

We have applied antisymmetrized molecular dynamics (AMD) with Gogny DIS force 
[19] to investigate the spectra of ^^Mg and ^^Ne. It will be shown that both nuclei have 
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the strongly deformed 3/2~ ground states with a neutron 3p2h configuration. Similar to 
= 19 system, different particle-hole configurations coexist within small excitation energy 

and their spin-parity are associated with the Nilsson-like orbits occupied by the last neutron. 

It is also shown that the one neutron separation energy of '^^Ne is rather small and ]?- shell 

neutron-halo structure is possibly formed. 

This paper is organized as follows. In the next section, the framework of AMD and 

the calculational procedure are explained. In the section III, the results are presented. 

The change of neutron particle-hole configurations as function of deformation is discussed. 

They are related to the spin-parity of the excitation spectra. The calculated spectra and 

transition probabilities are compared with the observations, and a theoretical assignment of 

the excitation spectra is suggested. The final section summarizes this work. 



II. THEORETICAL FRAMEWORK 



The applied theoretical method is the same as our previous work [11]. The deformed- 
basis AMD [20, 21] is combined with the generator coordinate method (GCM) to calculate 
the excitation spectra of = 21 system. The particle-hole configuration is investigated by 
the analysis of neutron single particle orbits. 



A. Wave function, Hamiltonian and variation 

The intrinsic wave function of the system with mass A is given by a Slater determinant 
of single particle wave packets, 

^int = A{ipi,(p2,-;VA}, ^i{r) = (j)i{r)xi^i, (1) 

where (fi is the ith single particle wave packet consisting of the spatial 0j, spin Xi and isospin 
parts. The local Gaussian located at Zj is employed as 0j, 

(pi{r) = exp|- Va(ra-^j=) 1, 

^ <r=x,y,z V J 

Xi = ttiXt + f^iXi, + lAP = 1, 

= proton or neutron. (2) 

Here Zj, ccj, f3i and i/g. are the variational parameters. The parity projected wave function, 
_ p^^.^^ is the variational wave function. 



The Gogny DIS force is employed as an effective nuclear force and the Coulomb force is 
approximated by a sum of seven Gaussians. 

In the present work, we have performed the variational calculation with the constraint on 
the matter quadrupole deformation /3. Variational parameters are optimized by the frictional 
cooling method so that the energy of the system is minimized for a given constraint on /3. 
The optimized wave function are denoted as It is noted that we do not put any 

constraint on the matter quadrupole deformation 7. Therefore 7 has the optimal value for 
each /3. 

B. Analysis of the single particle orbit 

To identify the particle-hole configuration, we analyze the neutron single particle orbital 
of <l>j„t(/3). We transform the single particle wave packet ipi to the orthonormalized basis, 
(fa — Ylit=i '^ia'fii- Here, Aq, and Qq, are the eigenvalues and eigenvectors of the overlap 
matrix B^j — {(pi\(fj). Using this basis, the Hartree-Fock single particle Hamiltonian, 

A 

7=1 




is defined. The eigenvalues and eigenvectors f^s of haf} give the single particle energies 
and the single particle orbits, 0s — Ylia=i fas^a- This procedure gives us information of the 
occupied neutron orbits. In this paper the particle-hole configuration is denoted as mpnh 
where m and n respectively denote number of neutron particles in p/-shell and holes in 
sd-shell in the spherical limit. For example, the normal configuration of A'^ = 21 system is 
denoted as IpOh. However, it is noted that this procedure only gives a rough estimation of 
the particle-hole configuration and should be regarded as a guide to understand the dominant 
particle-hole configuration of each states. The actual variational wave function is that 
is a sum of two Slater determinant. And we then perform the angular momentum projection 
and GCM calculation that superpose thousands of Slater determinants. Therefore our final 
wave function explained in the next subsection is beyond a simple single particle picture. 
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C. Angular momentum projection and GCM 



After the variational calculation, we project out an eigenstate of the total angular mo- 
mentum J, 

^mM-Pmk^HP) (4) 

Here P^k total angular momentum projector. The integrals over three Eulcr angles 

are evaluated by the numerical integration. The calculation is completed by performing 
GCM. The wave functions ^"j^xW) which have the same parity and angular momentum but 
have different K, deformation ^ and single particle configurations are superposed. Then the 
wave function of the system is expressed as 

= Cn^i^AP) + <*S.'(/30 + • • • , (5) 

where the quantum numbers except for the total angular momentum and the parity are 
represented by n. The coefficients c„, c^,... are determined by the Hill- Wheeler equation. 
The physical quantities discussed in the next section are basically calculated from the GCM 
wave function, Eq. (5). 

III. RESULTS AND DISCUSSIONS 

A. Energy curve and single particle configurations of ^^Mg 

Energy curves obtained by the variational calculation after the parity projection are 
presented in Fig. 1 (a) and (b) with circles. For a given value of (3, there is sometimes 
a local minimum solution together with the global minimum. In such case, both minima 
are plotted for the same value of /3. For any value of /3, ^^Mg prefered prolate deformation 
(7 > 30°) over oblate deformation (7 < 30°). The positive-parity states have two minima 
and they have neutron 2plh and Ap3h configurations in the ascending order of deformation. 
The negative-parity states have three minima with IpOh, 3p2h and 5pAh configurations. The 
particle-hole configurations of those minima are confirmed by the analysis of the neutron 
single particle orbits. Figure 1 (c) and (d) shows the single particle energies of the last 7 
neutrons as function of deformation. In the case of the negative parity, the last neutron 
occupies /7/2 orbit and other neutrons fill sd-shell at small deformation. Therefore, the 
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matter quadrupole deformation p matter quadrupole deformation p 



FIG. 1: (color online) (a) and (b): Energy curves of '^'^Mg before and after angular momentum 
projection for (a) positive- and (b) negative-parity states. Circles with dashed lines show the 
energy curve before the angular momentum projection. Colored lines show the energy curves after 
the angular momentum projection to J = 1/2, 3/2, 5/2 and 7/2. (c) and (d): Corresponding 
single particle energies of the last 7 neutrons. Filled (open) symbols show the orbits occupied 
by two (one) neutrons. Circles (boxes) show the orbits in which positive-parity (negative-parity) 
component is larger than 50 %. 

minimum around P = 0.2 in Fig. 1 (b) has a neutron IpOh configuration. As deformation 
becomes larger, neutron single particle energies show Nilsson model like behavior. Two 
orbits originate in fj/2 come down and splits into two orbits and go up. They are 
respectively denoted as [N,n„l„j,] = [3,3,0,1/2], [3,2,1,3/2], [2,0,2,3/2] and [2,0,0,1/2] 
in terms of Nilsson model. Around /3 = 0.25 the configuration changes. Two neutrons 
occupy [3,3,0,1/2] and one occupies [3,2,1,3/2] orbit. The [2,0,2,3/2] orbit is now empty. 
Thus the second minimum of the energy curve has a 3p2h configuration. Further increase 
of deformation pulls down the [3,2,1,1/2] orbit originates in ^3/2. The last neutron occupies 
it and two orbits originate in ^3/2 become empty. Therefore third minimum that appears 



6 



FIG. 2: Intrinsic density distributions at minima of energy curves, (a) and (b) show the minima 
of positive-parity states, (c) and (d) show Oplh and 3p2h minima of negative-parity states. 

as the local minimum above the 3p2h configuration has a 5pAh configuration. Its energy is 
rather higher than other two configurations. The positive-parity curve is also understood 
in the same way. It has two minima that have respectively 2plh and 4p3h configurations. 
The corresponding intrinsic matter density distributions are shown in Fig. 2. It shows an 
almost spherical shape of the normal configuration (IpO/i) and deformation becomes larger 
as number of neutrons in p/-shell increases. 

The energy curves and corresponding neutron single particle orbits show the qualitative 
agreement with the Nilsson model picture discussed in Ref. [22, 23]. The neutron single 
particle orbits split and change their order depending on deformation of the system. As 
deformation becomes larger, number of particle in the orbits originate in p/-shell increases. 
However, it is reminded that the analysis of the neutron single particle orbits gives only a 
rough estimation of the particle-hole configuration. Indeed, there are strong parity mixing of 
the neutron single particle orbits. For example, the orbit denoted as [3,2,1,1/2] in Fig. 1 (d) 



has the strongest parity mixing where the the positive-parity component amounts to about 
40%. The parity mixing becomes stronger as deformation larger. This may mean that the 
strongly deformed states have a complicated mixing of different particle-hole configurations 
and are beyond the Nilsson model picture. 

Adding to the Nilsson model like behavior, the present calculation has shown the fol- 
lowing two points. (1) ^^Mg is strongly deformed and quite soft against deformation. The 
intruder 3p2h configuration is already the lowest energy configuration at this stage and other 
particle-hole configurations also appear within very small excitation energy. (2) The last 
neutron walks around the [3,3,0,1/2], [2,0,2,3/2], [3,2,1,3/2], [2,0,0,1/2] and [3,2,1,1/2] orbits 
as deformation becomes larger. These neutron orbits will characterize the spin-parity of the 
corresponding excited states or rotational bands. 

Lines in Fig. 1 show the energy curves after the angular momentum projection. The 
energies of strongly deformed states are lowered by more than 5 MeV. The lowest energy state 
at this stage is the 3/2" state with the 3p2h configuration. More importantly, the property 
of the last neutron orbital is refiected to the excitation spectrum of each configuration. 
Deformation of IpOh configuration is the smallest among all minima and it generates low- 
lying 7/2~ and 3/2^ states that are respectively attributed to the almost spherical neutron 
/7/2 and P3/2 orbits. In the case of other minima, rotational spectra appear depending on 
the last neutron orbit. For example, the 3p2h configuration has the last neutron in the 
[3,2,1,3/2] orbit and it generates the K'^ — 3/2~ rotational spectrum. In the same way, the 
2pl/i, ApSh and 5p3h configurations respectively generate the K'^ — 3/2+, 1/2+ and l/2~ 
spectra that are associated with the asymptotic quantum number jz of the last neutron. 

B. Spectrum of ^^Mg and transitions 

Figure. 3 shows the low-lying spectrum of ^^Mg obtained by GCM. Tentative assignments 
given in Refs. [13, 15] are also shown. The calculated binding energy of ^^Mg is 251.1 MeV 
that underestimates the observed value by about 1 MeV. The one neutron separation energy 
is 1.9 MeV that is not far from the observed value (2.22 MeV). The ground state is the 3/2^ 
state and dominated by the 3p2h configuration that approximately amounts to 87%. Here 
the amount of the 3p2h configuration in the ground state is estimated by calculating the 
overlap between the ground state GCM wave function and the wave function of 3/2~ 
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FIG. 3: Calculated spectrum of ^^Mg compared with the tentative assignments suggested in Refs. 
[13, 15]. mpnh configurations shown below each bands mean the dominant particle configuration. 
Numbers show excitation energies in MeV. 



minimum 

($3/2 (/5 = 0.45)) in Fig. 1 (b), 



(6) 



The 3p2h configuration dominates the ground band K'" = 3/2" that has the first excited 
state 5/2i at 0.6 MeV. The normal configuration IpOh is located at higher excitation energy. 
The 7/22 and 3/22 states at 2.1 and 2.5 MeV are understood as the almost spherical fi/2 
and P3/2 single particle states. The 2plh and 3p2h configurations appear around 1 MeV. 
They respectively dominate the = 3/2"*" and 1/2"*" rotational bands. Similar to '^^Mg, 
many-particle and many-hole states with different deformation coexist within rather small 
excitation energy. 

The calculated B[E2) values are summarized in Table. I for several low-lying states. 
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FIG. 4: (color online) Calculated and observed [13, 14] /3 decay strength of ^^Na —?- ^^Mg (left) 
and ^^Mg ^^Al (right). 



Except for the IpOh configuration, the intra-band E2 transitions are large due to their 
deformation. Because of the mixing between different particle-hole configurations, there are 
non negligible inter-band E2 transitions between the K'^ = 3/2^*" and 1/2"'' bands. The P 
decay strength from ^^Na and to ^^Al are shown in Fig. 4. The wave functions of ^^Na and 
^^Al are calculated by using the same method as ^^Mg. The ground state of ^^Na is 3/2"'" 
and dominated by a neutron Ap2h configuration and mixed with a 2p0h configuration (10%) 
in the present calculation. Therefore the P decay strongly feeds the Ap3h configuration of 
^^Mg. The P decay to ^^Al mainly feeds the 3p3h configuration, since the ground state of 
33]\/[g is dominated by a 3p2h configuration. It feeds the 3/2~, 5/2"" and 1/2"" states around 
4.5 MeV strongly. 

Then we compare our results with observations and examine the spectrum of ^^Mg. 
Our result supports the 3/2^ assignment of the ground state suggested from the magnetic 
moment and two neutron knock-out reaction. Indeed the calculated 3/2^ state gives the 
closest value to the observed magnetic moment among the lowest energy states of particle- 
hole configurations (Table. II). The calculated first excited state (5/2^) at 0.6 MeV may 
correspond to the observed 0.48 MeV state, because the large Coulomb excitation cross 
section from the ground state to this state is reported and large quadrupole deformation 
Pc = 0.52 is extracted [24]. Experimentally, the first excited state was reported at 0.16 
MeV from the P decay measurement of '^'^Na. Since there was no direct P decay feeding 
of this state, a tentative assignment of 7/2~ was given. In our calculation, there is not 
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corresponding state with very small excitation energy. The observed 0.71 and 1.24 MeV 
states could be associated with the calculated 1/2+ (1.1 MeV) and 3/2^ (1.4 McV) states of 
K'^ = 1/2+ band from their log ft values. However, the present calculation does not explain 
the f3 decay of ^^Na to the ground state of ^^Mg nor that of ^^Mg to '^^Al. Especially, our 
calculation has no /3 decay feeding of low-lying states of ^^Al below 4 MeV, while there are 5 
transitions including the feeding of the ground state of ^^Al are reported. To explain those 
transitions, we need the positive-parity ground state of ^^Mg. 

An alternative assignment for the ground state of ^^Mg could be 1/2^^ with the ApSh 
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TABLE II: Calculated and observed [15] magnetic moments in the unit of /in. Theoretical values 
are shown for the lowest energy states of particle-hole configurations. 
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-0.91 


-0.7456 



configuration, since its magnetic moment has the same sign with and is not far from the 
observed value. Though the 1/2^ state is located at 1.1 MeV in the present calculation, a 
minor modification of the Gogny force parameter might make it the ground state. In this 
case, the observed /3 decay data of ^^Na could be compatible. The observed 0.71 and 1.24 
MeV states could be associated with the 3/22 and 5/22 states at 1.4 and 2.0 MeV from 
their log ft values. Their excitation energies measured from the 1/2^^ state are 0.4 and 1.0 
MeV. The observed 0.16 and 0.48 MeV states might be associated with the negative-parity 
states with IpOh or 3p2h configurations. However, even in this scenario, we cannot explain 
the /3 decay of ^^Mg. The l/2f state with ApSh configuration would strongly feed the states 
of ^^Al dominated by 4p4/i configuration that are located around 4.4 MeV in the present 
calculation. These transitions could be associated with the observed strong decays to 3.71, 
4.31 and 4.73 MeV states. The problem is that this assignment neither explains the feeding 
of the states below 4 MeV. There arc only the excited states with a 2p2h configuration 
and the ground state with the OpOh configurations below 4 MeV. Since there is always the 
mixing between different particle-hole configuration, the 1/2^ ground state of ^^Mg could 
feed these states, but it will never explain the strong feeding of the ground state of ^^Al that 
is dominated by OpOh configuration. 

Thus the present calculation partially explains some of observations but is incompatible 
with some in any interpretation. Further detailed study including the spectra of ^^Na and 
^^Al is required to settle down this problem, that will be done in near future. 

C. Spectrum of ^^Ne 

Since neutron number is the same, the energy curves and particle-hole configurations of 
^^Ne are quahtatively similar to ^^Mg. The difference is in the relative energies between 
particle- hole configurations. The calculated energy spectrum is shown in Fig. 5. The 
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FIG. 5: Calculated spectrum of ^^Ne. mpnh configurations shown below each bands mean the 
dominant particle configuration. Numbers show excitation energies in MeV. 

ground state is the 3/2~[ state dominated by the 3p2h configuration. Different from ^^Mg, 
the K'" = 3/2+ band {2plh) is located at smaller excitation energy and the K'" = 1/2+ band 
is at higher excitation energy. We conjecture that their excitation energies are sensitive to 
the neutron excess, since we have found a minor modification of the isospin dependence of 
Gogny force strongly affects their energies. 

Experimentally, small one neutron separation energy Sn = 0.29 ± 1.64 MeV is reported 
[25]. But neither the ground state spin-parity nor excitation spectrum are not well estab- 
lished. The assignment of ground state 3/2^ or 1/2+ is suggested from the Coulomb breakup 
reaction [18]. From these observations, the p-wave or s-wave neutron halo structure is sug- 
gested and theoretically investigated [26]. Therefore, we focus on the ground state property 
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of ^^Ne. The calculated one neutron separation energy is Sn — 250 keV that is fairly smaller 
than ^^Mg. Combined with the calculated ground spin-parity of 3/2~, ^^Ne is a promising 
candidate of a p-wave halo nucleus as suggested in Ref. [18]. The ground state of ^■'^Ne 
can be approximately regarded as the two body system consists of a ^°Ne core with 2p2h 
configuration and a neutron in [3,2,1,3/2] orbit. Since the ground band of ^°Ne is dominated 
by 2p2h configuration [27, 28], there should be the collective excitation of ^°Ne core in ^^Ne. 
The calculated matter root mean square radius of ^^Ne is 3.42 fm that is larger than ^°Ne 
(3.31 fm), ^^Ne (3.40 fm) and ^^Mg (3.34 fm). However, the neutron density distribution 
does not have a halo-like tail. This is due to the hmitation of the present calculation. Since 
the single particle wave function is limited to the Gaussian form, it cannot describe a halo 
tail. Improvement of the neutron wave function by combining the AMD with the resonating 
group method will be made in our next study. 

IV. SUMMARY 

In summary, we have studied the spectra of ^^Mg and ^^Ne. Both nuclei have the 
strongly deformed 3/2~ ground state with a 3p2h configuration and the coexistence of 
many different particle-hole configurations within small excitation energy. Their spectra 
are qualitatively understood in terms of the Nilsson model. The last neutron occupies 
[N,n,,l,J,] = [3,3,0,1/2], [2,0,2,3/2], [3,2,1,3/2], [2,0,0,1/2] and [3,2,1,1/2] depending 
on deformation of the system and the corresponding excited states or rotational bands are 
formed. The calculated spectrum of ^^Mg is examined based on the comparison with the 
observed magnetic moment and ^ decays. The 3/2~ assignment for the ground state by 
the present calculation agrees with the observed magnetic moment, but does not explain 
the f3 decays. An alternative assignment of 1/2"*" is also examined, but not compatible with 
the observed ^^Mg ^^Al (g.s.) decay. This requires further detailed study of the ^4 = 33 
isobars. ^^Ne also has the 3/2" ground state. The calculated small one neutron separation 
energy (250 keV) supports the p-wave neutron halo structure. 
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